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The Kinetic Sunyaev-Zel'dovich effect of the Milky Way Halo 
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We calculate the expected imprint of the ionized gas in the Milky- Way halo on the Cosmic Mi- 
crowave Background (CMB) through the kinetic Sunyaev-Zel'dovich (kSZ) effect. Unlike other 
Galactic foregrounds, the halo kSZ signature covers the full sky, generates anisotropies on large 
angular scales, is not accompanied by spectral distortions, and could therefore be confused with 
primordial CMB anisotropies. We construct theoretical models for various halo components, in- 
cluding smooth diffuse gas, filaments of cold inflowing gas and high velocity clouds. We find that 
the kSZ effect for all components is above the sensitivity of the Planck satellite, over a range of 
angular scales. However, the typical halo contribution is well below the cosmic variance noise in the 
primordial CMB power spectrum. High velocity clouds could dominate the halo contribution and 
better observational data is required to mask them out. We derive expected kSZ maps based on 
existing data from tracers of the halo gas distribution, such as 21cm maps of neutral hydrogen and 
H a maps of recombining gas. The cross-correlation of these maps with the WMAP5 data does not 
yield any statistically significant signal. 
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I. INTRODUCTION 

The scattering of the Cosmic Microwave Background 
(CMB) on free electrons with a bulk peculiar velocity rel- 
ative to the cosmic rest frame, shifts the observed CMB 
temperature and induces the so-called kinetic Sunyaev- 
Zel'dovich (kSZ) effect [kSZ,[IJ. Unlike its thermal coun- 
terpart which (for an isotropic distribution of electron ve- 
locities) alters the shape of the CMB spectrum, the kSZ 
effect shifts the spectrum while maintaining its black- 
body shape. This makes the kSZ signal on large angu- 
lar scales a foreground contamination that is indistin- 
guishable from the CMB anisotropies induced at the last 
scattering surface located at the redshift of cosmological 
recombination, z ~ 10 3 . 

The removal of the kSZ foreground requires knowledge 
of the distribution of free electrons and their velocity pat- 
tern. A local system of ionized gas that covers the entire 
sky and is capable of inducing kSZ anisotropies on large 
angular scales is the halo of our own Milky- Way (MW) 
galaxy. In this paper, we examine the kSZ signature im- 
printed by the ionized halo gas. Our primary goal is to 
find whether this unknown kSZ foreground is an impor- 
tant contaminant which limits our ability to extract the 
primordial anisotropies on the CMB sky. 

The MW galaxy can be modelled as a typical disk 
galaxy in the concordance Cold Dark Matter (ACDM) 
Universe. Theoretical studies of galaxy formation 
predict that in the most massive halos with > 1O 12 M , 
the infalling gas is shock-heated to the virial temperature 
of the halo, where it resides in hydrostatic equilibrium 
until cooling allows it to drift inwards and make stars in 
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a central galaxy. This configuration is obtained in hy- 
drodynamic simulations and observed in X-ray clusters. 
Lower-mass halos are too cold to emit hard X-rays, and 
it is challenging to observe their soft X-ray emission. The 
ionized gas in cluster halos causes a kSZ contamination of 
the CMB which is smaller than its thermal counterpart 
[{|. Dolag et al. [(J estimated the SZ effect of local large 
scale structures using constrained realization simulations 
and found the kSZ to be sub-dominant to the thermal 
SZ, both peaking at I ~ 1000. Suto et al. [7| suggested 
that the thermal effect of the local group might be de- 
tectable in CMB measurements [see, however, Ref. H, 
which claims that the assumed local group halo contra- 
dicts observations of other galaxy groups] . Other systems 
of ionized gas, including the intergalactic medium during 
the reionization era, also provides a negligible contribu- 
tion on large angular scales [for example, Ref. [9|, finds 
that the kSZ contribution from reionization peaks on the 
small angular scale of a few arcminutes] . 

Although the existence of gaseous halos around galax- 
ies is a robust theoretical prediction, its direct observa- 
tion has so far proven to be challenging [10]. Nicastro 
et al. [ll| and Rasmussen et al. |T2[ both find evidence 
of extended halo around the local group containing as 
much as 10 12 M Q and 10 9 — IO^Mq of gas respectively, 
by using OVII and OVIII absorption lines (and assum- 
ing a metallicity Z/Zp = 0.3-1). Other studies used the 
emission of OVII 13|, 14 1 , and the difference between the 
absorption line strength of the local group and the MW 
[lij to argue that the gas is more likely to be bound to 
the MW than the local group. This conclusion is con- 
sistent with Sembach et al. [15| . Yao et al. [16| which 
detect OVI lines but do not detect OVII line indicat- 
ing that the gas is bound to the shallower MW potential 
well. Hydrodynamic constraints from the interaction be- 
tween the Magellanic stream gas and the ambient dilute 
gas predict a non-negligible gaseous halo at a distance 
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of ~ 50kpc from the center of the MW [U El]. The 
existing results (see review in Ref. [11]) are sensitive to 
uncertain assumptions about the metalicity and the ex- 
istence of a multi-temperature or multi-component gas. 
We therefore follow a theoretical approach in modeling 
the gas in the MW halo and the local group. The lat- 
est models [H, [H| of the MW predict that it has a to- 
tal mass of M V i r i a i m 1 — 2 x 1O 12 M , a baryonic mass 
of M bar < 2 - 3 x lO n M , and a disk+bulge mass of 
Mdisk ~ 5 x 1O 1O M . This allows for a considerable frac- 
tion of the total baryonic mass to reside in a gaseous halo, 
partly in a hot phase at the virial temperature T « 10 6 K 
spread out to the virial radius i? V ir ~ 300kpc, and partly 
in a cooler halo component. 

Hydrodynamic simulations of galaxy formation predict 
that the hot MW halo is linked to its cosmological envi- 
ronment and should threfore be highly inhomogeneous. 
Typically, gas would be accreted into the halo along fil- 
aments which are oriented based on the MW's position 
within the cosmic web of dark matter and bary ons. Fol- 
lowing the results presented by Bregman 19] that the 
MW currently retains its separate halo, we expect fil- 
aments to be directed towards the MW's center rather 
than the Local group center. The gas in halos of galax- 
ies is not always heated by a virial shock near the virial 
radius. Birnboim & Dekel @, Q showed analytically that 
the post-shock gas is not always able to support the virial 
shock. If the halo mass is below some threshold, the post- 
shock gas is unstable to cooling, and the virial shock can- 
not propagate outwards. Rather, it remains close to the 
galactic disk. The threshold mass for metal-poor gas is 
~ 5 x lO n M and for gas at 0.3 of the solar metallicity 
it is ~ 1O 12 M . At high redshifts, halos above this mass 
threshold originate in rare density peaks, and such halos 
typically reside in the intersections of filaments within 
the cosmic web. This allows for denser filaments to avoid 
shock heating, even if the smooth, diffuse halo hosts a 
stable virial shock. At z > 2, halos above 1O 12 M are 
expected to contain cold flows within filaments embed- 
ded in a hot gaseous halo. This picture has been verified 
by different numerical simulations Keres et al. [4j, Ocvirk 
et al. HU and by H a observations [11]. The MW halo is 
just above the above mass threshold and it contains a hot 
halo. However, the filamentary structure is still expected 
to exist, as suggested by Keres et al. [24|, who used cos- 
mological Smooth Particle Hydrodynamics (SPH) simula- 
tions to show that overdense filaments of MW-type halos 
are more than an order of magnitude denser than the hot 
diffuse gas in which they are embedded, although 80-90% 
of the gas remains at the virial temperature. 

The hot gas in the MW halo has a cooling time that 
is shorter than the halo formation time interior to some 
radius [25] and its central core is expected to cool. De- 
tailed ID hydrodynamic simulations find, however, that 
the hot halo component does not cool significantly [11] 
for some period of time after the initial formation of the 
virial shock. As the post-shock gas becomes stable, the 
shock rapidly propagates from the central galaxy to the 



edge of the halo, halting the gas infall and, in some cases, 
putting it in an outward ballistic motion. This causes a 
complete shutdown of gas accretion which lasts for a few 
billion years. In the absence of a heating source, cool- 
ing would resume after this phase. The unimportance of 
"hot mode" accretion due to cooling of hot gas and set- 
tling to the galactic disk is also seen in 3D simulations of 
Keres et al. [24[. Mass limits on the galaxy indicate that 
even if gas had cooled, it did not settle into the galaxy, 
and would be orbiting the halo as cold clumps [but see 
Ref. H3, which demonstrated that the formation of cold 
clumps within the halo via a cooling fragmentation in- 
stability is highly unlikely] . In that case, the cold gas in 
the halo would be partially ionized [H, HI] and yield an 
even larger signature on the CMB anisotropies. 

A different local contribution to the kSZ effect can arise 
from high velocity clouds (HVCs). While the nature and 
formation scenarios of HVCs are still under debate [13] , 
we will adopt the inference (3(| that HVCs represent com- 
pact gas clouds with masses of 1O 6 -1O 8 M . The orbital 
velocities of HVCs are typically > 100 km s" 1 and they 
are typically found in the vicinity of the galactic plane 
(at \z\ < 10 kpc (although the latter may very well be 
an observational selection effect). Some of these clouds 
emit H a radiation [TtJ , indicating the existence of ionized 
hydrogen. The gas is ionized by sources in the Galactic 
disk, or by hydrodynamical heating (e.g. from the inter- 
action between the Magellanic stream and the hot halo 
gas [H, HH). The known distribution of free electrons 
in the Galactic disk yields a kSZ effect that is mostly 
confined to a narrow strip of the sky [H, HI], where 
other Galactic foregrounds (resulting from synchrotron 
and dust emission) dominate. 

In this work we calculate the kSZ effect from the ion- 
ized gas in the MW halo, which cover the entire sky. 
We proceed along two paths. First, we model in |JTT] 
three halo components: (i) a diffuse hot gas at the virial 
temperature; (ii) cold filaments; and (Hi) high velocity 
clouds, and then estimate the relative effects of these 
components on the CMB in £ )III1 Second, we adopt an 
empirical approach in £HV[ and cross-correlate various 
sky maps of halo gas tracers with the CMB data from 
WMAP5. Finally, we summarize our results and other 
possible sources of the kSZ effect in i|Vl 

II. KSZ SOURCES IN THE MILKY- WAY HALO 

Before examining the different kSZ components of the 
Galactic halo, let us obtain a rough estimate for their ex- 
pected magnitude. The kSZ shift in the CMB brightness 
temperature is given by [l|, 

AT f v R (D) 

— = - — - — a T n e {D)dD, (1) 

where vr is the radial peculiar velocity at a distance D 
along the line of sight (relative to the cosmic rest frame), 
n e {D) is the electron number density, c the speed of light 
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and ut is the Thomson cross section. A plausible col- 
umn density N e = 10 21 cnr 2 of halo electrons moving at 
velocity 600km s _1 (of order the velocity of the MW rela- 
tive to the CMB) yields a brightness temperature change 
of AT/T = 1.2 x 10~ 6 , which is detectable in current 
CMB measurements. 



A. The Milky- Way's invisible halo 

The baryonic mass fraction of the MW galaxy is ex- 
pected to be close to the average cosmic value, Qb/Q m s» 
0.167 ± 0.008 [ii], since the MW potential well is suffi- 
ciently deep (and the central black hole is of relatively low 
mass) for feedback to be inefficient in driving baryons out 
[35l | . Using highly efficient feedback, Dave |36| finds that 
65% of the gas remains in a MW sized halo. By subtract- 
ing the mass of stars and gas in the galaxy, 5-7 x 10 10 M Q 
[lO, [37j], from the total baryonic gas of ~ 1.5-3 x lO n M , 
we infer that the halo gas has a mass of > 10 11 M Q . 
Since the halo mass is larger than the critical mass for 
the formation of a virial shock [2j, |3j , a hydrostatic halo 
is expected to form. Detailed cosmological simulations 
predict that a sub-dominant fraction of the accreted 
baryons will be accreted in cold filaments which are pres- 
sure confined by the hot gaseous halo. This expectation 
is supported by Kaufmann et al. ^28| who find that under 
some initial conditions the halo gas fragments into cold 
and hot components, and that most of the gas remains 
in the halo for lOGyrs. We assume that at least half of 
the halo baryons, 5 x 10 10 M Q , reside in the hot halo. 

We assume that the total mass distribution is described 
by an NFW H profile with M mr = 1.5 x 10 12 M Q , cor- 
responding to the updated circular velocity of the MW 
[20l l2ll ]. Based on Refs. [13, H^|, we adopt a concen- 
tration parameter of c = 12. The virial radius for this 
mass is 296kpc. We consider two models for the halo gas 
distribution in hydrostatic equlibrium within the dark 
matter potential: (i) an isothermal (constant tempera- 
ture) profile ; and (ii) a generalized NFW density profile 

4C I with an inner core, in analogy with galaxy clusters 

41 1 and hydrodynamic simulations [42| , 
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where x = rc/ R vlr and p s is a normalization factor. In 
the isothermal models, the central density of the gas is 
set so that the total baryonic mass is encompassed within 
a radius where the baryonic overdensity is 30 (relative to 
the cosmic average) . For the core model, the gas density 
is normalized so as to give the total amount of assumed 
mass within i? V ir, and the temperature profile is fixed by 
the hydrostatic equilibrium equation. The resulting tem- 
perature profile is not far from isothermal, with an overall 
temperature variation by a factor of 2. Figure Q] shows 
the density and temperature profiles for the isothermal 
model (model "iso") with T iso = T vir = 8.1 x 10 5 .?r (solid 
red line), T iso = 1.2T vir = 9.7 x 10 5 K (model "1.2iso", 
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FIG. 1: Density (thick lines) and temperature (thin lines) of 
the virialized gas in the MW halo for three models: (i) an 
isothermal profile at the virial temperature T v i r (iso, red solid 
line); (ii) an isothermal profile at 1.2T v i r (1.2iso, long-dashed 
green line); and (Hi) NFW profile with a gaseous core (core, 
dashed blue line). The virial radius is 296kpc. 



long dashed green line) , and the core model (model core, 
dashed blue line). Although the profiles do not differ 
much from each other and they all include the same total 
baryonic mass within i? v ir, the electron column density 
of a radial line-of-sight from the center outwards changes 
significantly between N 2 i = (Ag/10 21 cm" 2 ) = 4.9, 1.7 
and 0.54 for models iso, 1.2iso and core, respectively. 

Kaufmann et al. [28| performed detailed SPH sim- 
ulations of MW-like galaxies based on previous work 
about cooling-induced fragmentation in galactic halos 
[29( 1 . They compared the halo cooling and galaxy accre- 
tion properties of two MW halo models: one with the gas 
in an NFW profile (their "low entropy" model), and the 
other with a core profile ("high entropy" model), much 
like our "core" model, with a density of n » 10~ 3 cm~ 3 
at a radius of 0.1i? V ir- They have found that while the 
first model cools the gas too fast, producing a disk with 
50% of the total baryonic mass, angular momentum defi- 
ciency, and a halo with a too high soft X-ray luminosity, 
the second model fragments via cooling instability into a 
configuration of cold clouds embedded within a hot dif- 
fuse halo which does not suffer from these problems. A 
fraction ~ 10% of the gas remains embedded in the hot 
halo as HVCs with a predicted covering factor of « 40%. 
The contribution of HVCs to the kSZ signature of the 
MW halo will be discussed in CfCl 

The rotation velocity of the hydrostatic MW halo is 
expected to be small, and so we assume that the halo is 
moving at a constant bulk velocity relative to the CMB 
0: Vgc-cmb = (552.2km s- 1 ,; = 266.5,6 = 29.1). 
We assume that the Sun is located at a distance of 8.5kpc 
from the galactic center [2(J HH, |3| • The off-center view- 
point creates a quadopolc, and higher spherical harmon- 
ics multipoles even for this simple configuration. 
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B. Filamentary infall into the Milky- Way halo 

In the ACDM cosmology, dark matter or gas are chan- 
neled into galactic halos through filaments which are in- 
terconnected to a cosmic web. The gas in the filaments is 
cold with a temperature in the range ~ 10 4 -10 5 K. As the 
cold gas enters the galactic halo, its subsequent thermo- 
dynamic evolution is determined by its metalicity and the 
gravitational potential well of its host halo. Previous an- 
alytic studies 0, [|| have found that in halos with a virial 
mass below 10 12 M Q and a metalicity of ~ 0.3-Z©, the gas 
is expected to fall to the center without shocking. For ha- 
los in a range just above this mass threshold at redshifts 
z > 2, a hot gaseous halo is expected to co-exist together 
with cold, unstable filaments. Detailed cosmological SPH 
simulations Q and Eulerian AMR simulations (22j in 3D, 
have confirmed these expectations. With the MW mass 
so close to the threshold mass, the filaments are still ex- 
pected to exist within the MW halo. The confinement 
of dense infalling filaments within the diffuse halo com- 
ponent was investigated in detail by Keres et al. [24l |. 
who have found that the filament gas is overdense by 1-2 
orders of magnitude relative to its hot environment and 
is infalling at ~ 200 km s _1 . The infalling gas stops at 
small radii in the vicinity of the galactic disk. Although 
the MW is expected to be fed by filaments, it is not pos- 
sible to predict their number, mass flux, and geometry. 
In this paper, we will assume that there are either one 
or three filaments, oriented radially towards the Galactic 
center, each having an infall speed of 200km s _1 , a coni- 
cal opening angle of 10° and an overdensity of 10 relative 
to the ambient hot halo. The direction of these filaments 
is chosen randomly; this particular choice has little effect 
on the kSZ power spectrum ( Alii Bp . 



C. High velocity clouds 

High velocity clouds (HVCs) are observed in 21cm HI 
data [30( as concentrations of gas which are moving at 
velocities |Vlsa| > 100km s _1 relative to the local stan- 
dard of rest (LSR) corresponding to the rotation of the 
Galaxy. The origin of the HVCs is still debated in the 
literature. Possible associations include the MW disk 
[see, e.g. Ref . |45|. and references therein], debris from the 
larg e or small Magellanic clouds [3l[ , extragalactic origin 
[27ll30j [but see also Ref.l46l which claims that the 21cm 
emission originates within a distance of 100 pc]. The dis- 
tance and ionization fraction of individual clouds (and 
therefore their size and mass) are unknown. Putman 
et al. [13 estimated distances of HVCs based on their 
H Q flux and models for the emission of ionizing radiation 
from the MW. They find, within the modelling and mea- 
surement uncertainties, that most HVCs are within a dis- 
tance of < 30kpc. Here we use these results to calculate 
the inferred electron column densities of HVCs by assum- 
ing, for simplicity, that they are spherical homogeneous 
clouds. For the angular diameter of HVCs on the sky, we 



adopt their typical observed diameter, 9 = 0.2rad [based 
on WHAM and LAB surveys; see Refs.SJQ- The total 
H Q photon flux / (in photons cm~ 2 s _1 sr) from a HVC 
at a distance D and angular diameter 9 is obtained from 
the recombination rate, 



/ = 



aBn 2 e VHVC 
4nD 2 



(3) 



where Vhvc = 47r /?>R z hvc is the volume of the cloud, 
with Rhvc = D x 9/2 being its physical radius. Here 
as = 1.17 x 10~ 13 sec -1 cm 3 is the case B recombination 
rate into H a photons . The electron column density 
through the center of the HVC is, 



N e = 2R HVC n e 



QnDf 



(4) 



In calculating the peculiar velocity of a HVC, we first 
transform the radial LSR velocity from Putman et al. 
[TtJ to the CMB frame of reference, based on the mea- 
surements of Kogut et al. [43| , 

Vhvc-cmb = Vhvc-lsr + Vlsr-gc + Vgc-cmb, (5) 

with Vhvc -LSR. = vrR where R is the unit vector in the 
HVC direction. The velocities have been transformed 
to Cartesian Galactic coordinates in which the x-axis is 
oriented towards the Galactic center and the z-axis is 
oriented perpendicular and out of the galactic plane. In 
these coordinates, Vlsr-cmb = (—37, —259, 269)km s _1 . 

We analyse 25 HVCs from Table 1 of Putman et al. pjj 
that include all the parameters needed. A more complete 
inventory of HVCs and their H Q measurements is neces- 
sary to produce foreground maps for future CMB exper- 
iments. The full contribution of these HVCs to the sky 
map will be shown in £11111 but we choose, for illustrative 
purpose, to consider the values from Putman et al. [l7| 
for a single HVC, "MSl-IIa" , which is a HVC in the Mag- 
ellanic stream. For this cloud vr-lsr — —120km s _1 
according to HI measurements and —124km s _1 accord- 
ing to H Q measurements. Converting to the CMB ref- 
erence frame, the peculiar radial velocity of this HVC 
is VhvC-Cmb X R = — 376km s _1 . For a distance 
D = 9.7kpc and an H Q flux of / = 407mR = 3.24 x 
10 4 photons cm -2 s^ 1 sr, we find N e = 1.7 x 10 20 cm -2 
and since the reported value of Njji is 1.1 x 10 20 cm -2 
we infer that this cloud is partly ionized. The mass of 
the cloud is 1.7 x lO 6 Af and its radius is 1 kpc. The 
kSZ effect on CMB photons passing through the center 
of the cloud is according to Ec^O, AT/T = 1.5 x 10~ 7 . 

The HVCs of Putman et al. [l7| were discovered by the 
proximity effect of the galactic ionizing radiation. This 



1 We use the table for case B, low density recombination rate from 
Osterbrock and Ferland [§4.2 of 49] for a gas temperature of 
10 4 K. 
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restricts their detections to the innermost 30kpc around 
the Galactic disk. Assuming that the origin on these 
HVCs is extra-galactic, there is no reason to assume that 
HVCs will only exist near the galactic disk. In addition, 
a comparison of the baryons in the galaxy and hot halo 
to the cosmic baryonic fraction of the total halo mass in- 
dicates that there are some missing baryons in the MW 
( QII A[ [28|, |29|). We therefore construct a synthetic dis- 
tribution of HVCs. We assume that they trace an isother- 
mal density profile, and include 2.5 x 10 10 M Q of mass, 
distributed as 250 HVCs, with a mass of 10 8 M Q each. 
Their direction is drawn randomly around the MWs cen- 
ter, and their distance from the MW is drawn from an 
isothermal density profile. The velocity has been drawn 
from a Gaussian distribution with a zero mean and an 
rms value, a = 100km s _1 . A constant angular size of 
6 = 0.2rad is assumed for all HVCs. The constant an- 
gular dimension is justified, for example, if the halo had 
an isothermal profile of gas and the HVCs at a different 
constant temperature (e.g. 10 4 K) were in hydrostatic 
equilibrium with the hot halo gas, and the displacement 
of the observer from the Galactic center is ignored. We 
have verified that a modest change in the radial distribu- 
tion and the angular size distribution of HVCs does not 
change their cumulative kSZ effect considerably. 



III. ESTIMATING THE KSZ EFFECT IN THE 
MW HALO 



We generate mock kSZ maps by first constructing a 3D 
model of the MW halo and then numerically integrating 
Eq. (TT|) over lines-of-sights in different directions. The 
3D density distribution of the halo electrons is split into 
a smooth spherically-symmetric component (see ^11 A[) . 
and "objects" which include cold filaments ( Hll Bl) and 
HVCs ( AH Cp . The various versions of the halo compo- 
nent are all static and centered around the Galactic cen- 
ter. Their column densities (from the Galactic center to 
the virial radius) change between 0.5 to4 x 10 21 cm -2 for 
the "core model" and the "iso model" , respectively. The 
filamets are conical objects, with an opening angle of 
10° and an overdensity of 10 with respect to the smooth 
halo. They are radially infalling at a constant speed of 
200km s _1 . Their central column densities are 10 times 
larger than that of the corresponding halo, but the ac- 
tual column density (integrated from the Sun) depends 
somewhat arbitrarily on their location on the sky (see 
QII B[) . The HVCs' velocity is either taken from obser- 
vations [ItJ or chosen randomly to fill the halo (typically 
\v r \ 100km s _1 ) and they have typical column densities 
of ~ 2 x 10 20 cm 2 ( jjTISJ) - The gas density and radial 
peculiar velocity at specific 3D positions are then inte- 
grated along each line-of-sight. We adopt a HEALPix 



[50t | 2 pixalization scheme with Nside= 256 which pro- 
vides sufficient resolution for the typical low multipolcs 
(I < 50) that we predict (see TO)) . The resulting kSZ 
sky maps are then decomposed into spherical harmon- 
ics, and the resulting power spectrum is compared to the 
CMB power spectrum, the cosmic variance noise, and 
the expected instrumental noise of the Planck satellite 3 . 
The analysis is performed using the Fortran 90 HEALPix 
software library. 



A. Sky maps 

We investigate a few halo profiles as described in Ta- 
ble |U Models 1-3 include only smooth profiles based on 
the halo models. Models 4 and 5 show the incremen- 
tal contribution of one filament at Galactic longitude I 
and latitude b values of (I = 259°, b = 30°) (model 4), 
and of two additional filaments at (I = 120°, b = 30°) 
and (I = 20°, b = -80°) (model 5) with respect to the 
galactic center. The locations of the filaments have been 
chosen arbitrarily for illustrative purposes. Model 6 adds 
to model 5 the HVC component according to Putman 
et al. ,17] ( mi CI) that has a sub-dominant effect rela- 
tive to the halo and filament contributions. Model 7 de- 
scribed a synthetic HVC population of 250 HVCs with 
total mass of 2.5 x 10 10 M o as described in mi CI Model 
8 is the maximal effect model achieved by combining a 
hydrostatic iso profile, 3 filaments and the synthetic dis- 
tribution of filaments. The maximal effect of model 8 is 
AT/T - 8 x 10~ 6 , which exceeds the values of 2[iK re- 
ported by Waelkens et al. [32 |. Hajian et al. [33| for the 
Galactic disk. The square of the signal-to-noise ratio is 
defined as, 

= Tl h-TT) £ f™&< ^ 

\ l max l min\ i -} ^ -t V ^ 

with l min = 4, l max = 20 and N^ MAP5 being the to- 
tal error estimate of the unbinned WMAP5 data 4 . We 
show the associated S/N values in the last column of Ta- 
ble HI The models with synthetic HVCs in Table U have 
(S/N) ~ 0.2, making them sub-dominant relative to cos- 
mic variance, although not by much. If the real clouds 
have a higher column density of electrons than assumed 
here (e.g. owing to a non-spherical geometry), then their 
contribution could become noticeable. 

Figure [5] shows the kSZ sky maps for a smooth halo 
model, a halo model with filaments, and a halo model 
with filaments and HVCs (models 3, 5 and 7, respec- 
tively) with the monopole and dipole components owing 



2 http://healpix.jpl.nasa.gov 

3 http:/ /www.rssd.esa.int/index.php?projcct=planck 

4 taken from http://lambda.gsfc.nasa.gov/data/map/dr3/dcp/ 
wmap_tt_spectrum_5yr_v3pl .txt 
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TABLE I: the kSZ effects of various halo models. In the HVCs 
column '-','+','*' corresponds to no HVCs, HVCs according 
to Putman et al. [TtJ, and synthetic HVCs, respectively. 



# 


halo 
model 


number of 
filaments 


HVCs 


min(^-) 

[io-*] 


mai(-y-) 

[io- d ] 


S/N 


1 


ISO 







-3.0 


3.5 


.030 
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to our motion removed. The dipole resulting from the 
relative velocity of the halo electrons and the CMB is 
shown. 



B. kSZ Power spectra 

The power spectrum is recovered from the sky maps 
by decomposing the sky maps into spherical harmonics, 
[AT/T] (6, <j>) = J2 hm a^°Y lm (9,cf>), and averaging over 

m to obtain C* al ° = 2IT1 £m \ a \m°\ 2 for each *• Fig- 
ure 0] compares the (binned) power spectra of the 5-year 
WMAP data [H to that of models 3,5,6 and 7. For 
reference, the figure also shows the statistical noise aris- 
ing from cosmic variance. The figure also presents the 
expected instrumental noise of the Planck satellite, 

Nr p = (m,m) = 9Lw4exp{-^ + }, (7) 

with ri; the noise of each o; m , 6f w hm = 10' the full width 
at half maximum of the beam and <jt = bfiK being the 
detector's sensitivity 5 , and the instrumental noise as re- 
ported in the WMAP 5 data. 

We find that the effect of the ionized MW halo is larger 
than Planck's instrumental noise for I < 4. The power 
spectrum scales roughly linearly with the number of fila- 
ments. The filament contribution peaks at I ~ 15 and is 
larger than Planck's instrumental noise for I < 20. The 
observed HVCs contribution to the power spectrum is 
sub-dominant relative to that of the filaments. The syn- 
thetic HVCs contribute to the power spectrum at I ~ 20 
and are larger than Planck's instrumental noise at I < 40. 
A more complete observational census of HVCs with H a 
measurements (beyond the 25 HVCs known) will increase 
the signal linearly in the number of HVCs and improve 
these constraints. All of the kSZ power spectra lie well 



5 The values quoted here are reasonable for the 143GHz HFI 
instrument, http:/ /www.rssd.esa.int/SA/PLANCK/docs/goal- 
perf-sum-feb-2004.pdf 
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FIG. 2: Fractional changes in the CMB brightness temper- 
ature, AT/T, from the kSZ effect of various models for the 
MW halo. All three halo models include a smooth component 
of the "core" density profile, described in O Al and Fig. [T] 
Upper panel: the smooth core model (model 3). Middle panel: 
core model with three filaments (model 5). Lower panel: core 
model with three filaments and synthetic HVCs (model 7). 



below the cosmic variance noise, making them statisti- 
cally insignificant for the determination of cosmological 
parameters. All of the models discussed here are below 
the instrumental noise of WMAP5, and not expected to 
be observed. 

The kSZ contributions from filaments and HVCs are 
localized and should not in general resemble a Gaus- 
sian random field. In an attempt to quantify the non- 
Gaussianity of the kSZ signal we generated mock CMB 
maps using the power spectrum of the WMAP5 [HI and 
superimposed the results of model 7 onto these maps. 
We then calculated the brightness fluctuations on dif- 
ferent angular scales around the HEALPix cell centers 
and compared a histogram of these to a corresponding 
histogram of the mock CMB map without the filament 
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FIG. 3: The one point Gaussianity histogram of random 
maps, and random maps with model 7 (core, 3 filaments and 
250 HVC) superimposed on it. The histogram is generated by 
polling the average value of 8T on disks with diameter 0.15rad. 
The shaded area is the la scatter of 20 such realizations. The 
blue lines correspond to the non-Gaussian contribution when 
model 7 is multiplied by 1, 10 and 30. The normalization is 
arbitrary. 



and HVC contributions. For various choices of the an- 
glular scale and the binning intervals, we have found 
that the non-Gaussian kSZ signature of HVCs and fil- 
aments is undetectable. If the number of the HVCs is 
increased by an order of magnitude, deviations from a 
Gaussian probability distribution are still virtually non- 
detectable. Figure [3] compares the non-Gaussianity of 
model 7, when the kSZ signal is multiplied by factors 
of 1, 10, and 30. We quantify the strength of the non- 
Gaussian signal by measuring the "power in the wings" 
which we define as the ratio between the integral of the 
histogram beyond some threshold fluctuation amplitude 
for the non-Gaussian and Gaussian distributions. The 
cutoff chosen here is ^4 x 10~ 5 K (thin vertical line). 
The power in the wings is: 1.0025,1.066 and 1.439 for 
the factors of 1, 10, and 30 respectively. We conclude 
that if the signal were an order of magnitude larger (a 
factor of 10 in the column density of the ionized gas or the 
velocity), then the non-Gaussianity would have been de- 
tectable. A more rigorous three-point correlation (Jnl) 
test goes beyond the scope of this paper. 



IV. CROSS CORRELATIONS OF HI AND H a 
MAPS WITH THE CMB 

In the previous sections we utilized physically moti- 
vated models for the halo gas including filaments and 
HVCs. In this section we follow an empirical approach 
and cross correlate various observational maps that might 
trace the distribution of free electrons in the MW halo. 
The procedure used here follows the analysis of Land and 
Slosar [5l| for the Leiden/Argentine/Bonn (LAB) Galac- 



1Q- 2 - «..<'»v'"» 



10- 




1 (multipole) 

FIG. 4: The CMB WMAP5 power spectrum (red curve with 
error bars) in comparison to the power spectrum of the core 
halo profile (model 3, long-dashed green line), core with 3 fil- 
aments (model 5, dotted blue line) and core with 3 filaments 
and HVCs (model 6, purple dashed line) and core + 3 fila- 
ments + synthetic HVCs (model 7, cyan dot-dashed). The 
thin solid lines delineate the cosmic variance uncertainties 
(higher line) and Planck's instrumental noise (lower). The 
instrumental noise of the WMAP5 data is denoted by the 
thin, dashed line. Features in the shaded area below the in- 
strumental noise line will not be observable with the Planck 
satellite. 



tic HI Survey [521 ]. with the appropriate density/ velocity 
weighting as discussed below. We start by describing the 
cross correlation of a general sky map (that may cover 
only part of the sky) with the CMB in §IVAI We then 
analyse various specific maps and weighting schemes in 
3TVCJ1 



A. Cross-correlation formalism 

We adopt a procedure for cross-correlating a sky map 
with the CMB data that resembles the one used by Land 
and Slosar (5lj | , with the important difference that we use 
a weighting scheme tailored to enhance the kSZ signal. 
Our procedure is as follows: 

i. Interpolate the sky map into a HEALPix format 
with Nside= 512. HI maps originating from the 
LAB survey are in fits image format, and so we 
performed a 2D bilinear interpolation to the coor- 
dinates of the center of each HEALPix cell. For H Q 
WHAM maps [47J , the original data is presented as 
scatter plots and we use qshed2d algorithm [53| to 
interpolate it into HEALPix cell positions. 

ii. Remove the monopole and dipole of the maps by 
converting the map to (I, m) space and eliminating 
all the power in I < 1, and then converting back to 
real space. 

iii. Combine the mask of the sky map with the 
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WMAP5 data mask KQ85 [H to produce the "to- 
tal mask" , and convolve the total mask with both 
WMAP data (ILC and V2 map were used) and the 
sky map. 

iv. Calculate the cross-correlation of the WMAP5 map 
with the sky map under consideration using the 
total mask. The auto/cross correlation of maps a 
and b is given by: 

c r 6 = (aTi)E a ^ ( g ) 

and the normalized cross correlation is defined as: 

s~ta,b 

C?> h = , 1 ■ (9) 

A perfect correlation yields C ; a ' b — 1, a null corre- 
lation yields 0, and a perfect anti-correlation yields 
a value of —1 for this parameter. 

v. Since the real maps are expected to be non- 
Gaussian and we expect the nontrivial masks to 
introduce some correlations between different mul- 
tipoles, we evaluate the statistical significance of 
the result by creating random mock CMB maps 
from the unmasked CMB power-spectrum, and cal- 
culating the masked, normalized cross-correlation 
10, 000 times. The mean and variance of the real- 
izations for each / are then compared to the actual 
cross-correlation signal. 



B. Observed Sky maps 

LAB survey: The kSZ weight should depend linearly on 
the electron column density and the velocity of the ion- 
ized gas (Eq. [T]). We use data on the neutral h ydr ogen 
distribution from the 21cm LAB survey [52l 1551 l56l l57l]. 
which provides a full sky map of the column density at 
each velocity bin relative to the LSR. The data is tabu- 
lated in bins of 0.5° in angle and a velocity resolution of 
1.3 km s — 1 . The data covers the full sky, and so we use 
the mask of KQ85 which leaves 82% of the sky unmasked. 

The conversion of the measured HI column density JVm 
to the column density of free electrons N e requires some 
model assumptions. We first consider the case where 
the gas is optically thin to ionizing radiation and is dis- 
tributed in a spherical shell of a constant width D around 
the galaxy. Ionizing radiation, either extragalactic or 
galactic in origin, is impinging on the gas and creating 
free electrons. In this steady state balance be- 

tween the recombination and ionization rates imply, 

n\D oc n m D = N m , (10) 

with n e , rifji and Nhi denoting the electron density, the 
HI density and the HI column density, respectively. The 



LAB survey maps Ahi as a function of velocity and so 
the sky map LABI is constructed as 

= - / N H1 (l, b, vf^vcMBdv, (11) 

J v 

with ucmb = v + Vlsr-cmb x f, Vlsr-cmb from Kogut 
et al. [43J and f the radius vector to the point of interest 
(l,b). 

As an alternative case, we assume that the gas is op- 
tically thick to ionizing radiation and is collisionally ion- 
ized at some constant temperature. In this case, 

n e D(xnffiD = N H i, (12) 

and 

= - / N H i(l,b,v)v CM Bdv. (13) 

J v 

Outliers in HI: In addition to cross correlating the full 
HI LAB data, we also correlated the CMB with a derived 
map of HVCs and outliers produced by Westmeier [58| 
by subtracting an HI halo model [5!| from the data, and 
masking all the areas with no extra MW component. The 
outlier map masks all but 12% of the sky and after con- 
volving this mask with the CMB mask, the total mask 
leaves 10% of the sky unmasked. The data is expressed 
in terms of the average velocity and total column density 
and in analogy with Eqs. (fTTj) and (TT3")) we define models 
OL1 and OL2 as, 

S? b L1 = -<*(!, b)v CM B, (14) 

and 

S°b 2 = -N H i(l,b)v C MB. (15) 

H a data: The HI data, while being extremely accu- 
rate and detailed, is an indirect tracer of free electrons. 
A more direct tracer is H Q radiation which result from 
recombinations of electrons and protons. The Wiscon- 
sin H- Alpha Mapper [WHAM, provides an H Q map 
of the sky. Above the galactic plane, the WHAM sur- 
vey has been shown to correlate with HVCs [17] . The 
WHAM survey covers 73% of the sky and the total mask 
amounts to 61% of the sky. Assuming, again, some typ- 
ical width for the H Q emitting layer, D, and using Eq. 
Q , the column density of ionized gas should scale as the 
square root of the H Q flux, f 1 ^ 2 , and so we adopt the 
weighting scheme, 

STi HAM = - f f(l,b,v) 1/2 v CMB dv. (16) 

The WHAM data is available as velocity dependent emis- 
sion measure at scattered locations, and was interpolated 
to the HEALPix grid following Renka [53] . 
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FIG. 5: Normalized cross correlation as a function of spherical 
harmonic multipoles. We show the full LAB data with LABI 
weighting (Eq. 1111 top panel), with LAB2 weighting (Eq. 
[H 2nd panel), HVC map with LABI (3rd panel), HVC with 
LAB2 (4th panel), and WHAM with weighting according to 
Eq. (|16[l . Solid line: cross correlation between each map and 
WMAP5 [3||- Shaded area: 1-cr scatter levels from 10,000 
mock CMB realizations with the WMAP5 power spectrum. 
Data has been binned to groups of 5 consecutive multipoles in 
order to reduce the noise level. The X 2 of the cross correlation 
and the fraction of mock realizations with \ 2 larger than the 
real one are marked on each panel. 




(5lT | have performed a similar analysis on the LAB data, 
with sky maps of column density of integrals over all 
or part of velocity space (rather than our expected kSZ 
weightings), and have found no significant cross correla- 
tion. We recover these results here and find for our LABI 
and LAB2 sky maps that there is no apparent cross cor- 
relation with the WMAP5 data. A null cross correlation 
is also consistent with the OL2 and WHAM sky maps. 

Sky map OL1 shows an anti correlation larger than 
la over most multipoles between 50 < I < 200. The 
statistical significance of this anti correlation was tested 
by defining 



(17) 



with l m i n = 50, l m ax = 200 and 07 being the standard 
deviation of 10 4 realization. The \ 2 °f the actual data 
was then compared to the distribution of \ 2 for 10 re- 
alizations. 6% of the mock realizations had a x 2 larger 
than the observed sky map, indicating that such an out- 
come is not an extremely rare occurrence. Additionally, 
we performed various tests to eliminate possible biases, 
including rotation of the OL map by an arbitrary angle 
(which recovers the null result and eliminates some effects 
of the mask), and checking different frequency bands of 
WMAP data (and finding that the ILC and V bands yield 
almost similar, statistically significant anti-correlation). 
The anti correlation is probably due to correlations of 
different multipoles introduced by the nontrivial mask 
which covers 89% of the sky, while a true kSZ detection 
should have yielded a positive correlation. We also note 
that the OL sky maps were derived from the Westmeier 
[HI maps that were produced by subtracting an HI halo 
model from the LAB data. An over subtraction could 
have produced an anti correlation of the type that was 
found. Finally, the LAB, WHAM and outlier maps have 
been searched (non-systematically) for features that re- 
semble some of the predicted halo features in iJTT] and no 
clear correlation was found. 



DISCUSSION AND CONCLUSIONS 



C. Cross correlation results 

Figure [5] shows the normalized cross correlation (Eq. 

of the sky maps defined in ^IVBl The powers here 
have been binned to groups of 5 Is in order to reduce 
the noise level. A null cross correlation should have a 
zero correlation among different multipoles so this aver- 
aging should not change the expected (zero) mean of the 
cross correlation signal. The shaded areas on the plots 
mark the averaged la regions around averaged means of 
10, 000 cross correlations for each sky map with a mock 
CMB map of the same power spectrum. Land and Slosar 



We examined the kSZ effect from the MW halo by 
first constructing theoretical models for the expected sig- 
nal, and then adopting an empirical strategy of cross- 
correlation observed maps of halo gas tracers with the 
WMAP5 data. Our theoretical models include smooth 
hot gas, filaments of cold inflowing g ctS, clS well as HVCs. 
The cross-correlation was analysed for the observed sky 
maps of HI and H Q emission. 

Figure [J] shows that the kSZ signatures from the 
smooth halo and the filaments declines sharply for mul- 
tipoles above I = 4 and 20, respectively, and is unlikely 
to be detectable considering the orders of magnitude gap 
between the signal and the CMB power spectrum. The 
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contribution from HVCs peaks at I ~ 20 with an am- 
plitude which is ~ 3 orders of magnitude lower than the 
CMB power-spectrum and ~ 2 orders of magnitude lower 
than the statistical noise of WMAP5. However, for mul- 
tipoles I < 20-40, the total kSZ signal from the Galac- 
tic halo is expected to be above the instrumental noise 
of the Planck satellite. Galactic foregrounds are being 
routinely masked or cleaned in the analysis procedure, 
often based on their unique dependence on photon fre- 
quency or Galactic latitude. However, the kSZ effect of 
the MW halo covers the entire sky and has no frequency 
dependence relative to the blackbody spectrum, and so 
its removal cannot be accomplished in the standard lin- 
ear combination (ILC) maps. The halo effect is larger 
by an order of magnitude than the kSZ signature of the 
MW disk [33, HH , and is not confined to a narrow strip 
around the galactic plane - which is often masked out 
in a cosmological analysis of the CMB anisotropics. Our 
theoretical prediction for the HVC contribution to the 
kSZ effect is based on H Q observations of these HVCs 
[T?! together with a number of model assumptions. In 
particular, we assume the HVCs to be spherical objects 
with constant angular size, optically thin, and in ioniza- 
tion/recombination equilibrium with the galactic radia- 
tion field. Better understanding of the nature, location 
and geometry of these HVCs is required for a more robust 
assessment of their kSZ effect. The modest sample of 25 
HVCs with known distances [17J is incomplete, and since 
the kSZ signal depends linearly on the number of HVCs, 
a more complete data base is necessary. We estimated 
the effects of 250 such HVCs by randomly picking them 
at larger radii. The HVCs are expected to leave a distinct 
non-Gaussian signature on the CMB sky. Although our 
preliminary tests have found this signature to be statis- 
tically insiginificant, the future discovery of more HVCs 
with known distances might add a non-negligible non- 
Gaussian contamination to the primordial CMB signal. 

The kSZ and thermal SZ effects caused by large scale 
structures in our local environment has been tested by 
Dolag et al. [6J] using constrained realization simulations. 
They have found a kSZ signal comparable to the one pre- 
sented in this paper, but at larger multipoles. The ther- 
mal component they predict is even larger. The local 
group was proposed as a possible source for the thermal 
Sunyaev-Zel'dovich (tSZ) effect 0. The early proposals 
for the dipole and quadropole tSZ amplitudes of the lo- 
cal group are unlikely to be real based on measurements 
of halos of other galaxy groups The kSZ signature 
of the local group would lead to a compact spot in the 
direction of the Andromeda galaxy (i.e. the direction of 
the center of mass of the local group) , with a diameter of 
~ 20°. Unfortunately, the detection of this spot is highly 
challenging. 



Our empirically-calibrated sky maps of the kSZ signal 
are based on the HI data [52j , the HI HVCs and outliers 
distribution (58[, and the H Q emission maps [47]. We 
examined two algorithms for converting HI maps to elec- 
tron column densities: (i) gas which is optically-thin to 
ionizing radiation, and (ii) gas with a constant ionization 
fraction. We have found no statistically significant cross- 
correlation between the HI and H Q data and the CMB 
data. 

We conclude that the kSZ signal from HVCs imprints 
pattens on the CMB sky which are detectable but are sta- 
tistically insignificant in the cosmological context. The 
kSZ effect could potentially be used to study the content 
of the Galactic halo. A detection of non-Gaussianity with 
the spectral features predicted here (I < 20) would corre- 
spond to objects that are either very big or within the lo- 
cal group. Fine tuning of the kSZ effect could be invoked 
to explain some non-Gaussian features that are seen in 
the WMAP data. The cold spot [6(| can be explained 
by a dense HVC cloud receding from us at a relatively 
high speed. The details of such a cloud and the statisti- 
cal probability that such a cloud exists is left for further 
investigation. 

While cross correlation between various observable sky 
maps and CMB data is routinely done, we derive here a 
prescription to estimate the specific contribution from 
the kSZ effect, that incorporates the column density and 
velocity in a self consistent way. If a sky map traces 
the distribution of ionized gas, but its line of sight veloc- 
ities are not highly correlated, a cross correlation between 
that map and the CMB will fail to detect the kSZ feature 
unless the velocities are takes into account in a manner 
described in flIVl Reversely, a situation in which a sky 
map convolved according to fllVI prescriptions that pro- 
duces a correlation while the naive luminosity or column 
density map fails to yield one would indicate that kSZ is 
observed at high probability. Finally, a more complete 
future compilation of HVC data could help in creating 
masks that remove this foreground contamination as well 
as constrain the unknown physical properties of HVCs. 
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